################################################################################
# Libraries and Import
################################################################################

library(readstata13)
library(ggplot2)

### This file shows how the data set is prepared based on the raw Qualtrics response data.

dat <- read.dta13("data/China_Data_FINAL.dta")

################################################################################
# Data Cleaning
################################################################################

# incompletes: 309 dropped when downloading from Qualtrics

# early terminates 116 due to age and foreign residence
dat <- dat[dat$q27 != "", ]

# code vignette group
dat$id = 1:nrow(dat)
dat$vignette_group = NA
dat$vignette_group[dat$q52 != ""] <- 1
dat$vignette_group[dat$q62 != ""] <- 2
dat$vignette_group[dat$q72 != ""] <- 3
dat$vignette_group[dat$q82 != ""] <- 4

dat <- dat[!is.na(dat$vignette_group), ]
  
### Demographics

dat$age_group <- factor(dat$q22, 
       levels = c("20-24", "25-29", "30-34", "35-39", "40-44",
                  "45-49", "50-54", "55-59", "60及60以上"), 
       labels = c(1,2,3,4,5,6,7,8,9))

dat$male <- 0
dat$male[dat$q23=="男性"] <- 1

dat$prov <- as.factor(dat$q24)

dat$temp <- dat$q25
dat$temp[dat$temp == "其他"] <- "农村"
dat$current_resi <- as.factor(dat$temp)
dat$temp <- NULL

dat$temp <- dat$q26
dat$temp[dat$temp == "其他"] <- "农村"
dat$past_resi <- as.factor(dat$temp)
dat$temp <- NULL

# education; merge primary/secondary with high school
dat$edu <- dat$q27
dat$edu[dat$edu == "小学或以下"] <- "高中"
dat$edu[dat$edu == "初中"] <- "高中"
dat$edu <- factor(dat$edu,
                  levels = c("高中", "大专", "大学本科", "研究生或以上学历"), 
                  labels = c(1,2,3,4))
  
# regime_insider
dat$regime_insider <- 0
dat$regime_insider[dat$q28 == "国有企事业单位"| dat$q28 == "党政机关"] <- 1

# party member
dat$party <- 0
dat$party[dat$q29 == "共产党员"|dat$q29 == "共青团员"] <- 1

# income level
dat$income <- factor(dat$q210, 
                        levels = c("少于10,000", "10,001~30,000元", "30,001-60,000 元", "60,001-90,000 元", 
                                   "90,001-120,000 元", "120,001-200,000元",
                                   "高于200,000元"), 
                        labels = c(1,2,3,4,5,6,7))

# Japan exporter/importer
dat$Japan_export <- 0
dat$Japan_export[dat$q213 == "我所在的公司主要对日本出口商品"] <- 1 

dat$Japan_import <- 0
dat$Japan_import[dat$q213 == "我所在的公司主要从日本进口商品"] <- 1 

# dependency with Japan
dat$China_dep <- 0
dat$China_dep[dat$q37 == "中国更为依赖从日本进口的商品"|dat$q37 == "中国更为依赖向日本出口商品"] <- 1 

dat$Mutual <- 0
dat$Mutual[dat$q37 == "中日两国在贸易关系上互相依赖"] <- 1 

dat$Japan_dep <- 0
dat$Japan_dep[dat$q37 == "日本更为依赖从中国进口的商品"|dat$q37 == "日本更为依赖向中国出口商品"] <- 1

dat$Japan_influence <- 0
dat$Japan_influence[dat$q38 == "我会受到比较负面的经济影响"|dat$q38 == "我会受到非常负面的经济影响"] <- 1

### Political Attitudes
dat$perception_fp <- 0
dat$perception_fp <- factor(dat$q32, 
                            levels = c("非常强硬", "比较强硬", "介于强硬和温和之间",
                                       "比较温和", "非常温和"), 
                            labels = c(5,4,3,2,1))

dat$temp <- dat$q33
dat$temp[dat$temp == "不存在显著的影响"] <- "正面和负面的影响互相抵消"
dat$international_sys <- factor(dat$temp, levels = c("非常正面的影响", "比较正面的影响",
                                           "正面和负面的影响互相抵消", "比较负面的影响",
                                           "非常负面的影响"),
                                labels = c(5,4,3,2,1))
dat$temp <- NULL

dat$territory <- factor(dat$q34, levels = c("强烈赞同", "赞同",
                                            "既不赞同也不反对", "反对",
                                            "强烈反对"),
                        labels = c(5,4,3,2,1))

dat$China_citizen <- factor(dat$q35, levels = c("强烈赞同", "赞同",
                                             "既不赞同也不反对", "反对",
                                             "强烈反对"),
                            labels = c(5,4,3,2,1))

dat$critique <- factor(dat$q36, levels = c("强烈赞同", "赞同",
                                           "既不赞同也不反对", "反对",
                                           "强烈反对"),
                            labels = c(5,4,3,2,1))

# The last one 
dat$q310[dat$q310 == "不清楚/拒绝回答"] <- "否"
dat$fight <- factor(dat$q310, levels = c("是 ", "否"),
                    labels = c(1, 0))

################################################################################
# Format vignette data for analysis
################################################################################
temp <- data.frame(matrix("", nrow = nrow(dat), ncol = 14), stringsAsFactors=FALSE)

temp[dat$vignette_group == 1, 1:14] <- dat[dat$vignette_group == 1, which(colnames(dat) == "q52"):which(colnames(dat) == "q58_8")]
temp[dat$vignette_group == 2, ] <- dat[dat$vignette_group == 2, which(colnames(dat) == "q62"):which(colnames(dat) == "q68_8")]
temp[dat$vignette_group == 3, ] <- dat[dat$vignette_group == 3, which(colnames(dat) == "q72"):which(colnames(dat) == "q78_8")]
temp[dat$vignette_group == 4, ] <- dat[dat$vignette_group == 4, which(colnames(dat) == "q82"):which(colnames(dat) == "q88_8")]

colnames(temp) <- c("D1", "D3", "D4", "D5", "D6", "D7", 
                    "D8_a", "D8_b", "D8_c", "D8_d", "D8_e",
                    "D8_f", "D8_g", "D8_h")

dat_vignette <- cbind(dat, temp)

dat_vignette$D1 <- factor(dat_vignette$D1, levels = c("强烈赞同", "赞同",
                                           "既不赞同也不反对", "反对",
                                           "强烈反对"),
                       labels = c(5,4,3,2,1))

dat_vignette$D3 <- factor(dat_vignette$D3, levels = c("强烈赞同", "赞同",
                                          "既不赞同也不反对", "反对",
                                          "强烈反对"),
                       labels = c(5,4,3,2,1))

dat_vignette$D4 <- factor(dat_vignette$D4, levels = c("强烈赞同", "赞同",
                                    "既不赞同也不反对", "反对",
                                    "强烈反对"),
                 labels = c(5,4,3,2,1))

dat_vignette$D7 <- factor(dat_vignette$D7, levels = c("国际社会对中国政府的决策抱有正面的态度",
                                    "国际社会对中国政府的决策抱有中性的态度",
                                    "国际社会对中国政府的决策抱有负面的态度"),
                 labels = c(1,2,3))

dat_vignette$D8_a <- factor(dat_vignette$D8_a, levels = c("比较可能",
                                    "比较不可能"),
                 labels = c(1,0))

dat_vignette$D8_b <- factor(dat_vignette$D8_b, levels = c("比较可能",
                                        "比较不可能"),
                   labels = c(1,0))

dat_vignette$D8_c <- factor(dat_vignette$D8_c, levels = c("比较可能",
                                        "比较不可能"),
                   labels = c(1,0))

dat_vignette$D8_d <- factor(dat_vignette$D8_d, levels = c("比较可能",
                                        "比较不可能"),
                   labels = c(1,0))

dat_vignette$D8_e <- factor(dat_vignette$D8_e, levels = c("比较可能",
                                        "比较不可能"),
                   labels = c(1,0))

dat_vignette$D8_f <- factor(dat_vignette$D8_f, levels = c("比较可能",
                                        "比较不可能"),
                   labels = c(1,0))

dat_vignette$D8_g <- factor(dat_vignette$D8_g, levels = c("比较可能",
                                        "比较不可能"),
                   labels = c(1,0))

dat_vignette$D8_h <- factor(dat_vignette$D8_h, levels = c("比较可能",
                                        "比较不可能"),
                   labels = c(1,0))

dat_vignette[, c(which(colnames(dat) == "q42"):which(colnames(dat) == "q88_8"), 
                 which(colnames(dat) == "f11"):which(colnames(dat) == "f327"))] <- NULL

dat_vignette_sub <- dat_vignette[dat_vignette$q31== "其他", ]


################################################################################
# Format conjoint data for analysis
################################################################################
dat$conjoint1 <- as.numeric(dat$q42 == "政策１")
dat$conjoint1_a <- as.numeric(dat$q43_1)
dat$conjoint1_b <- as.numeric(dat$q43_2)

dat$conjoint2 <- as.numeric(dat$q45 == "政策１")
dat$conjoint2_a <- as.numeric(dat$q46_1)
dat$conjoint2_b <- as.numeric(dat$q46_2)

dat$conjoint3 <- as.numeric(dat$q48 == "政策１")
dat$conjoint3_a <- as.numeric(dat$q49_1)
dat$conjoint3_b <- as.numeric(dat$q49_2)


dat_conjoint <- data.frame(id = c(rep(dat$id[dat$f11 != ""], 2), 
                                   rep(dat$id[dat$f21 != ""], 2),
                                   rep(dat$id[dat$f31 != ""], 2)),
                            profile = c(rep(1, sum(dat$f11 != "")), 
                                        rep(2, sum(dat$f11 != "")),
                                        rep(1, sum(dat$f21 != "")), 
                                        rep(2, sum(dat$f21 != "")),
                                        rep(1, sum(dat$f31 != "")),
                                        rep(2, sum(dat$f31 != ""))))

dat_conjoint$chosen <- c(dat$conjoint1[dat$f11 != ""], as.numeric(dat$conjoint1[dat$f11 != ""] == 0),
                          dat$conjoint2[dat$f21 != ""], as.numeric(dat$conjoint2[dat$f21 != ""] == 0),
                          dat$conjoint3[dat$f31 != ""], as.numeric(dat$conjoint3[dat$f31 != ""] == 0))

dat_conjoint$score <- c(dat$conjoint1_a[dat$f11 != ""], dat$conjoint1_b[dat$f11 != ""], 
                         dat$conjoint2_a[dat$f21 != ""], dat$conjoint2_b[dat$f21 != ""],
                         dat$conjoint3_a[dat$f31 != ""], dat$conjoint3_b[dat$f31 != ""])

dat$C1a靖国神社 <- ""
dat$C1b靖国神社 <- ""
dat$C2a靖国神社 <- ""
dat$C2b靖国神社 <- ""
dat$C3a靖国神社 <- ""
dat$C3b靖国神社 <- ""

dat$C1a靖国神社[dat$f11=="靖国神社"] <- dat$f111[dat$f11=="靖国神社"]
dat$C1b靖国神社[dat$f11=="靖国神社"] <- dat$f121[dat$f11=="靖国神社"]
dat$C1a靖国神社[dat$f12=="靖国神社"] <- dat$f112[dat$f12=="靖国神社"]
dat$C1b靖国神社[dat$f12=="靖国神社"] <- dat$f122[dat$f12=="靖国神社"]
dat$C1a靖国神社[dat$f13=="靖国神社"] <- dat$f113[dat$f13=="靖国神社"]
dat$C1b靖国神社[dat$f13=="靖国神社"] <- dat$f123[dat$f13=="靖国神社"]
dat$C1a靖国神社[dat$f14=="靖国神社"] <- dat$f114[dat$f14=="靖国神社"]
dat$C1b靖国神社[dat$f14=="靖国神社"] <- dat$f124[dat$f14=="靖国神社"]
dat$C1a靖国神社[dat$f15=="靖国神社"] <- dat$f115[dat$f15=="靖国神社"]
dat$C1b靖国神社[dat$f15=="靖国神社"] <- dat$f125[dat$f15=="靖国神社"]
dat$C1a靖国神社[dat$f16=="靖国神社"] <- dat$f116[dat$f16=="靖国神社"]
dat$C1b靖国神社[dat$f16=="靖国神社"] <- dat$f126[dat$f16=="靖国神社"]
dat$C1a靖国神社[dat$f17=="靖国神社"] <- dat$f117[dat$f17=="靖国神社"]
dat$C1b靖国神社[dat$f17=="靖国神社"] <- dat$f127[dat$f17=="靖国神社"]

dat$C2a靖国神社[dat$f21=="靖国神社"] <- dat$f211[dat$f21=="靖国神社"]
dat$C2b靖国神社[dat$f21=="靖国神社"] <- dat$f221[dat$f21=="靖国神社"]
dat$C2a靖国神社[dat$f22=="靖国神社"] <- dat$f212[dat$f22=="靖国神社"]
dat$C2b靖国神社[dat$f22=="靖国神社"] <- dat$f222[dat$f22=="靖国神社"]
dat$C2a靖国神社[dat$f23=="靖国神社"] <- dat$f213[dat$f23=="靖国神社"]
dat$C2b靖国神社[dat$f23=="靖国神社"] <- dat$f223[dat$f23=="靖国神社"]
dat$C2a靖国神社[dat$f24=="靖国神社"] <- dat$f214[dat$f24=="靖国神社"]
dat$C2b靖国神社[dat$f24=="靖国神社"] <- dat$f224[dat$f24=="靖国神社"]
dat$C2a靖国神社[dat$f25=="靖国神社"] <- dat$f215[dat$f25=="靖国神社"]
dat$C2b靖国神社[dat$f25=="靖国神社"] <- dat$f225[dat$f25=="靖国神社"]
dat$C2a靖国神社[dat$f26=="靖国神社"] <- dat$f216[dat$f26=="靖国神社"]
dat$C2b靖国神社[dat$f26=="靖国神社"] <- dat$f226[dat$f26=="靖国神社"]
dat$C2a靖国神社[dat$f27=="靖国神社"] <- dat$f217[dat$f27=="靖国神社"]
dat$C2b靖国神社[dat$f27=="靖国神社"] <- dat$f227[dat$f27=="靖国神社"]

dat$C3a靖国神社[dat$f31=="靖国神社"] <- dat$f311[dat$f31=="靖国神社"]
dat$C3b靖国神社[dat$f31=="靖国神社"] <- dat$f321[dat$f31=="靖国神社"]
dat$C3a靖国神社[dat$f32=="靖国神社"] <- dat$f312[dat$f32=="靖国神社"]
dat$C3b靖国神社[dat$f32=="靖国神社"] <- dat$f322[dat$f32=="靖国神社"]
dat$C3a靖国神社[dat$f33=="靖国神社"] <- dat$f313[dat$f33=="靖国神社"]
dat$C3b靖国神社[dat$f33=="靖国神社"] <- dat$f323[dat$f33=="靖国神社"]
dat$C3a靖国神社[dat$f34=="靖国神社"] <- dat$f314[dat$f34=="靖国神社"]
dat$C3b靖国神社[dat$f34=="靖国神社"] <- dat$f324[dat$f34=="靖国神社"]
dat$C3a靖国神社[dat$f35=="靖国神社"] <- dat$f315[dat$f35=="靖国神社"]
dat$C3b靖国神社[dat$f35=="靖国神社"] <- dat$f325[dat$f35=="靖国神社"]
dat$C3a靖国神社[dat$f36=="靖国神社"] <- dat$f316[dat$f36=="靖国神社"]
dat$C3b靖国神社[dat$f36=="靖国神社"] <- dat$f326[dat$f36=="靖国神社"]
dat$C3a靖国神社[dat$f37=="靖国神社"] <- dat$f317[dat$f37=="靖国神社"]
dat$C3b靖国神社[dat$f37=="靖国神社"] <- dat$f327[dat$f37=="靖国神社"]

dat_conjoint$靖国神社 <- c(dat$C1a靖国神社[dat$f11 != ""], dat$C1b靖国神社[dat$f11 != ""],
                        dat$C2a靖国神社[dat$f21 != ""], dat$C2b靖国神社[dat$f21 != ""],
                        dat$C3a靖国神社[dat$f31 != ""], dat$C3b靖国神社[dat$f31 != ""])


dat$C1a经济合作 <- ""
dat$C1b经济合作 <- ""
dat$C2a经济合作 <- ""
dat$C2b经济合作 <- ""
dat$C3a经济合作 <- ""
dat$C3b经济合作 <- ""

dat$C1a经济合作[dat$f11=="经济合作"] <- dat$f111[dat$f11=="经济合作"]
dat$C1b经济合作[dat$f11=="经济合作"] <- dat$f121[dat$f11=="经济合作"]
dat$C1a经济合作[dat$f12=="经济合作"] <- dat$f112[dat$f12=="经济合作"]
dat$C1b经济合作[dat$f12=="经济合作"] <- dat$f122[dat$f12=="经济合作"]
dat$C1a经济合作[dat$f13=="经济合作"] <- dat$f113[dat$f13=="经济合作"]
dat$C1b经济合作[dat$f13=="经济合作"] <- dat$f123[dat$f13=="经济合作"]
dat$C1a经济合作[dat$f14=="经济合作"] <- dat$f114[dat$f14=="经济合作"]
dat$C1b经济合作[dat$f14=="经济合作"] <- dat$f124[dat$f14=="经济合作"]
dat$C1a经济合作[dat$f15=="经济合作"] <- dat$f115[dat$f15=="经济合作"]
dat$C1b经济合作[dat$f15=="经济合作"] <- dat$f125[dat$f15=="经济合作"]
dat$C1a经济合作[dat$f16=="经济合作"] <- dat$f116[dat$f16=="经济合作"]
dat$C1b经济合作[dat$f16=="经济合作"] <- dat$f126[dat$f16=="经济合作"]
dat$C1a经济合作[dat$f17=="经济合作"] <- dat$f117[dat$f17=="经济合作"]
dat$C1b经济合作[dat$f17=="经济合作"] <- dat$f127[dat$f17=="经济合作"]

dat$C2a经济合作[dat$f21=="经济合作"] <- dat$f211[dat$f21=="经济合作"]
dat$C2b经济合作[dat$f21=="经济合作"] <- dat$f221[dat$f21=="经济合作"]
dat$C2a经济合作[dat$f22=="经济合作"] <- dat$f212[dat$f22=="经济合作"]
dat$C2b经济合作[dat$f22=="经济合作"] <- dat$f222[dat$f22=="经济合作"]
dat$C2a经济合作[dat$f23=="经济合作"] <- dat$f213[dat$f23=="经济合作"]
dat$C2b经济合作[dat$f23=="经济合作"] <- dat$f223[dat$f23=="经济合作"]
dat$C2a经济合作[dat$f24=="经济合作"] <- dat$f214[dat$f24=="经济合作"]
dat$C2b经济合作[dat$f24=="经济合作"] <- dat$f224[dat$f24=="经济合作"]
dat$C2a经济合作[dat$f25=="经济合作"] <- dat$f215[dat$f25=="经济合作"]
dat$C2b经济合作[dat$f25=="经济合作"] <- dat$f225[dat$f25=="经济合作"]
dat$C2a经济合作[dat$f26=="经济合作"] <- dat$f216[dat$f26=="经济合作"]
dat$C2b经济合作[dat$f26=="经济合作"] <- dat$f226[dat$f26=="经济合作"]
dat$C2a经济合作[dat$f27=="经济合作"] <- dat$f217[dat$f27=="经济合作"]
dat$C2b经济合作[dat$f27=="经济合作"] <- dat$f227[dat$f27=="经济合作"]

dat$C3a经济合作[dat$f31=="经济合作"] <- dat$f311[dat$f31=="经济合作"]
dat$C3b经济合作[dat$f31=="经济合作"] <- dat$f321[dat$f31=="经济合作"]
dat$C3a经济合作[dat$f32=="经济合作"] <- dat$f312[dat$f32=="经济合作"]
dat$C3b经济合作[dat$f32=="经济合作"] <- dat$f322[dat$f32=="经济合作"]
dat$C3a经济合作[dat$f33=="经济合作"] <- dat$f313[dat$f33=="经济合作"]
dat$C3b经济合作[dat$f33=="经济合作"] <- dat$f323[dat$f33=="经济合作"]
dat$C3a经济合作[dat$f34=="经济合作"] <- dat$f314[dat$f34=="经济合作"]
dat$C3b经济合作[dat$f34=="经济合作"] <- dat$f324[dat$f34=="经济合作"]
dat$C3a经济合作[dat$f35=="经济合作"] <- dat$f315[dat$f35=="经济合作"]
dat$C3b经济合作[dat$f35=="经济合作"] <- dat$f325[dat$f35=="经济合作"]
dat$C3a经济合作[dat$f36=="经济合作"] <- dat$f316[dat$f36=="经济合作"]
dat$C3b经济合作[dat$f36=="经济合作"] <- dat$f326[dat$f36=="经济合作"]
dat$C3a经济合作[dat$f37=="经济合作"] <- dat$f317[dat$f37=="经济合作"]
dat$C3b经济合作[dat$f37=="经济合作"] <- dat$f327[dat$f37=="经济合作"]

dat_conjoint$经济合作 <- c(dat$C1a经济合作[dat$f11 != ""], dat$C1b经济合作[dat$f11 != ""],
                        dat$C2a经济合作[dat$f21 != ""], dat$C2b经济合作[dat$f21 != ""],
                        dat$C3a经济合作[dat$f31 != ""], dat$C3b经济合作[dat$f31 != ""])


dat$C1a中国对日本的关税水平 <- ""
dat$C1b中国对日本的关税水平 <- ""
dat$C2a中国对日本的关税水平 <- ""
dat$C2b中国对日本的关税水平 <- ""
dat$C3a中国对日本的关税水平 <- ""
dat$C3b中国对日本的关税水平 <- ""

dat$C1a中国对日本的关税水平[dat$f11=="中国对日本的关税水平"] <- dat$f111[dat$f11=="中国对日本的关税水平"]
dat$C1b中国对日本的关税水平[dat$f11=="中国对日本的关税水平"] <- dat$f121[dat$f11=="中国对日本的关税水平"]
dat$C1a中国对日本的关税水平[dat$f12=="中国对日本的关税水平"] <- dat$f112[dat$f12=="中国对日本的关税水平"]
dat$C1b中国对日本的关税水平[dat$f12=="中国对日本的关税水平"] <- dat$f122[dat$f12=="中国对日本的关税水平"]
dat$C1a中国对日本的关税水平[dat$f13=="中国对日本的关税水平"] <- dat$f113[dat$f13=="中国对日本的关税水平"]
dat$C1b中国对日本的关税水平[dat$f13=="中国对日本的关税水平"] <- dat$f123[dat$f13=="中国对日本的关税水平"]
dat$C1a中国对日本的关税水平[dat$f14=="中国对日本的关税水平"] <- dat$f114[dat$f14=="中国对日本的关税水平"]
dat$C1b中国对日本的关税水平[dat$f14=="中国对日本的关税水平"] <- dat$f124[dat$f14=="中国对日本的关税水平"]
dat$C1a中国对日本的关税水平[dat$f15=="中国对日本的关税水平"] <- dat$f115[dat$f15=="中国对日本的关税水平"]
dat$C1b中国对日本的关税水平[dat$f15=="中国对日本的关税水平"] <- dat$f125[dat$f15=="中国对日本的关税水平"]
dat$C1a中国对日本的关税水平[dat$f16=="中国对日本的关税水平"] <- dat$f116[dat$f16=="中国对日本的关税水平"]
dat$C1b中国对日本的关税水平[dat$f16=="中国对日本的关税水平"] <- dat$f126[dat$f16=="中国对日本的关税水平"]
dat$C1a中国对日本的关税水平[dat$f17=="中国对日本的关税水平"] <- dat$f117[dat$f17=="中国对日本的关税水平"]
dat$C1b中国对日本的关税水平[dat$f17=="中国对日本的关税水平"] <- dat$f127[dat$f17=="中国对日本的关税水平"]

dat$C2a中国对日本的关税水平[dat$f21=="中国对日本的关税水平"] <- dat$f211[dat$f21=="中国对日本的关税水平"]
dat$C2b中国对日本的关税水平[dat$f21=="中国对日本的关税水平"] <- dat$f221[dat$f21=="中国对日本的关税水平"]
dat$C2a中国对日本的关税水平[dat$f22=="中国对日本的关税水平"] <- dat$f212[dat$f22=="中国对日本的关税水平"]
dat$C2b中国对日本的关税水平[dat$f22=="中国对日本的关税水平"] <- dat$f222[dat$f22=="中国对日本的关税水平"]
dat$C2a中国对日本的关税水平[dat$f23=="中国对日本的关税水平"] <- dat$f213[dat$f23=="中国对日本的关税水平"]
dat$C2b中国对日本的关税水平[dat$f23=="中国对日本的关税水平"] <- dat$f223[dat$f23=="中国对日本的关税水平"]
dat$C2a中国对日本的关税水平[dat$f24=="中国对日本的关税水平"] <- dat$f214[dat$f24=="中国对日本的关税水平"]
dat$C2b中国对日本的关税水平[dat$f24=="中国对日本的关税水平"] <- dat$f224[dat$f24=="中国对日本的关税水平"]
dat$C2a中国对日本的关税水平[dat$f25=="中国对日本的关税水平"] <- dat$f215[dat$f25=="中国对日本的关税水平"]
dat$C2b中国对日本的关税水平[dat$f25=="中国对日本的关税水平"] <- dat$f225[dat$f25=="中国对日本的关税水平"]
dat$C2a中国对日本的关税水平[dat$f26=="中国对日本的关税水平"] <- dat$f216[dat$f26=="中国对日本的关税水平"]
dat$C2b中国对日本的关税水平[dat$f26=="中国对日本的关税水平"] <- dat$f226[dat$f26=="中国对日本的关税水平"]
dat$C2a中国对日本的关税水平[dat$f27=="中国对日本的关税水平"] <- dat$f217[dat$f27=="中国对日本的关税水平"]
dat$C2b中国对日本的关税水平[dat$f27=="中国对日本的关税水平"] <- dat$f227[dat$f27=="中国对日本的关税水平"]

dat$C3a中国对日本的关税水平[dat$f31=="中国对日本的关税水平"] <- dat$f311[dat$f31=="中国对日本的关税水平"]
dat$C3b中国对日本的关税水平[dat$f31=="中国对日本的关税水平"] <- dat$f321[dat$f31=="中国对日本的关税水平"]
dat$C3a中国对日本的关税水平[dat$f32=="中国对日本的关税水平"] <- dat$f312[dat$f32=="中国对日本的关税水平"]
dat$C3b中国对日本的关税水平[dat$f32=="中国对日本的关税水平"] <- dat$f322[dat$f32=="中国对日本的关税水平"]
dat$C3a中国对日本的关税水平[dat$f33=="中国对日本的关税水平"] <- dat$f313[dat$f33=="中国对日本的关税水平"]
dat$C3b中国对日本的关税水平[dat$f33=="中国对日本的关税水平"] <- dat$f323[dat$f33=="中国对日本的关税水平"]
dat$C3a中国对日本的关税水平[dat$f34=="中国对日本的关税水平"] <- dat$f314[dat$f34=="中国对日本的关税水平"]
dat$C3b中国对日本的关税水平[dat$f34=="中国对日本的关税水平"] <- dat$f324[dat$f34=="中国对日本的关税水平"]
dat$C3a中国对日本的关税水平[dat$f35=="中国对日本的关税水平"] <- dat$f315[dat$f35=="中国对日本的关税水平"]
dat$C3b中国对日本的关税水平[dat$f35=="中国对日本的关税水平"] <- dat$f325[dat$f35=="中国对日本的关税水平"]
dat$C3a中国对日本的关税水平[dat$f36=="中国对日本的关税水平"] <- dat$f316[dat$f36=="中国对日本的关税水平"]
dat$C3b中国对日本的关税水平[dat$f36=="中国对日本的关税水平"] <- dat$f326[dat$f36=="中国对日本的关税水平"]
dat$C3a中国对日本的关税水平[dat$f37=="中国对日本的关税水平"] <- dat$f317[dat$f37=="中国对日本的关税水平"]
dat$C3b中国对日本的关税水平[dat$f37=="中国对日本的关税水平"] <- dat$f327[dat$f37=="中国对日本的关税水平"]

dat_conjoint$中国对日本的关税水平 <- c(dat$C1a中国对日本的关税水平[dat$f11 != ""], dat$C1b中国对日本的关税水平[dat$f11 != ""],
                        dat$C2a中国对日本的关税水平[dat$f21 != ""], dat$C2b中国对日本的关税水平[dat$f21 != ""],
                        dat$C3a中国对日本的关税水平[dat$f31 != ""], dat$C3b中国对日本的关税水平[dat$f31 != ""])




dat$C1a日本对中国的关税水平 <- ""
dat$C1b日本对中国的关税水平 <- ""
dat$C2a日本对中国的关税水平 <- ""
dat$C2b日本对中国的关税水平 <- ""
dat$C3a日本对中国的关税水平 <- ""
dat$C3b日本对中国的关税水平 <- ""

dat$C1a日本对中国的关税水平[dat$f11=="日本对中国的关税水平"] <- dat$f111[dat$f11=="日本对中国的关税水平"]
dat$C1b日本对中国的关税水平[dat$f11=="日本对中国的关税水平"] <- dat$f121[dat$f11=="日本对中国的关税水平"]
dat$C1a日本对中国的关税水平[dat$f12=="日本对中国的关税水平"] <- dat$f112[dat$f12=="日本对中国的关税水平"]
dat$C1b日本对中国的关税水平[dat$f12=="日本对中国的关税水平"] <- dat$f122[dat$f12=="日本对中国的关税水平"]
dat$C1a日本对中国的关税水平[dat$f13=="日本对中国的关税水平"] <- dat$f113[dat$f13=="日本对中国的关税水平"]
dat$C1b日本对中国的关税水平[dat$f13=="日本对中国的关税水平"] <- dat$f123[dat$f13=="日本对中国的关税水平"]
dat$C1a日本对中国的关税水平[dat$f14=="日本对中国的关税水平"] <- dat$f114[dat$f14=="日本对中国的关税水平"]
dat$C1b日本对中国的关税水平[dat$f14=="日本对中国的关税水平"] <- dat$f124[dat$f14=="日本对中国的关税水平"]
dat$C1a日本对中国的关税水平[dat$f15=="日本对中国的关税水平"] <- dat$f115[dat$f15=="日本对中国的关税水平"]
dat$C1b日本对中国的关税水平[dat$f15=="日本对中国的关税水平"] <- dat$f125[dat$f15=="日本对中国的关税水平"]
dat$C1a日本对中国的关税水平[dat$f16=="日本对中国的关税水平"] <- dat$f116[dat$f16=="日本对中国的关税水平"]
dat$C1b日本对中国的关税水平[dat$f16=="日本对中国的关税水平"] <- dat$f126[dat$f16=="日本对中国的关税水平"]
dat$C1a日本对中国的关税水平[dat$f17=="日本对中国的关税水平"] <- dat$f117[dat$f17=="日本对中国的关税水平"]
dat$C1b日本对中国的关税水平[dat$f17=="日本对中国的关税水平"] <- dat$f127[dat$f17=="日本对中国的关税水平"]

dat$C2a日本对中国的关税水平[dat$f21=="日本对中国的关税水平"] <- dat$f211[dat$f21=="日本对中国的关税水平"]
dat$C2b日本对中国的关税水平[dat$f21=="日本对中国的关税水平"] <- dat$f221[dat$f21=="日本对中国的关税水平"]
dat$C2a日本对中国的关税水平[dat$f22=="日本对中国的关税水平"] <- dat$f212[dat$f22=="日本对中国的关税水平"]
dat$C2b日本对中国的关税水平[dat$f22=="日本对中国的关税水平"] <- dat$f222[dat$f22=="日本对中国的关税水平"]
dat$C2a日本对中国的关税水平[dat$f23=="日本对中国的关税水平"] <- dat$f213[dat$f23=="日本对中国的关税水平"]
dat$C2b日本对中国的关税水平[dat$f23=="日本对中国的关税水平"] <- dat$f223[dat$f23=="日本对中国的关税水平"]
dat$C2a日本对中国的关税水平[dat$f24=="日本对中国的关税水平"] <- dat$f214[dat$f24=="日本对中国的关税水平"]
dat$C2b日本对中国的关税水平[dat$f24=="日本对中国的关税水平"] <- dat$f224[dat$f24=="日本对中国的关税水平"]
dat$C2a日本对中国的关税水平[dat$f25=="日本对中国的关税水平"] <- dat$f215[dat$f25=="日本对中国的关税水平"]
dat$C2b日本对中国的关税水平[dat$f25=="日本对中国的关税水平"] <- dat$f225[dat$f25=="日本对中国的关税水平"]
dat$C2a日本对中国的关税水平[dat$f26=="日本对中国的关税水平"] <- dat$f216[dat$f26=="日本对中国的关税水平"]
dat$C2b日本对中国的关税水平[dat$f26=="日本对中国的关税水平"] <- dat$f226[dat$f26=="日本对中国的关税水平"]
dat$C2a日本对中国的关税水平[dat$f27=="日本对中国的关税水平"] <- dat$f217[dat$f27=="日本对中国的关税水平"]
dat$C2b日本对中国的关税水平[dat$f27=="日本对中国的关税水平"] <- dat$f227[dat$f27=="日本对中国的关税水平"]

dat$C3a日本对中国的关税水平[dat$f31=="日本对中国的关税水平"] <- dat$f311[dat$f31=="日本对中国的关税水平"]
dat$C3b日本对中国的关税水平[dat$f31=="日本对中国的关税水平"] <- dat$f321[dat$f31=="日本对中国的关税水平"]
dat$C3a日本对中国的关税水平[dat$f32=="日本对中国的关税水平"] <- dat$f312[dat$f32=="日本对中国的关税水平"]
dat$C3b日本对中国的关税水平[dat$f32=="日本对中国的关税水平"] <- dat$f322[dat$f32=="日本对中国的关税水平"]
dat$C3a日本对中国的关税水平[dat$f33=="日本对中国的关税水平"] <- dat$f313[dat$f33=="日本对中国的关税水平"]
dat$C3b日本对中国的关税水平[dat$f33=="日本对中国的关税水平"] <- dat$f323[dat$f33=="日本对中国的关税水平"]
dat$C3a日本对中国的关税水平[dat$f34=="日本对中国的关税水平"] <- dat$f314[dat$f34=="日本对中国的关税水平"]
dat$C3b日本对中国的关税水平[dat$f34=="日本对中国的关税水平"] <- dat$f324[dat$f34=="日本对中国的关税水平"]
dat$C3a日本对中国的关税水平[dat$f35=="日本对中国的关税水平"] <- dat$f315[dat$f35=="日本对中国的关税水平"]
dat$C3b日本对中国的关税水平[dat$f35=="日本对中国的关税水平"] <- dat$f325[dat$f35=="日本对中国的关税水平"]
dat$C3a日本对中国的关税水平[dat$f36=="日本对中国的关税水平"] <- dat$f316[dat$f36=="日本对中国的关税水平"]
dat$C3b日本对中国的关税水平[dat$f36=="日本对中国的关税水平"] <- dat$f326[dat$f36=="日本对中国的关税水平"]
dat$C3a日本对中国的关税水平[dat$f37=="日本对中国的关税水平"] <- dat$f317[dat$f37=="日本对中国的关税水平"]
dat$C3b日本对中国的关税水平[dat$f37=="日本对中国的关税水平"] <- dat$f327[dat$f37=="日本对中国的关税水平"]

dat_conjoint$日本对中国的关税水平 <- c(dat$C1a日本对中国的关税水平[dat$f11 != ""], dat$C1b日本对中国的关税水平[dat$f11 != ""],
                        dat$C2a日本对中国的关税水平[dat$f21 != ""], dat$C2b日本对中国的关税水平[dat$f21 != ""],
                        dat$C3a日本对中国的关税水平[dat$f31 != ""], dat$C3b日本对中国的关税水平[dat$f31 != ""])


dat$C1a日本修改宪法 <- ""
dat$C1b日本修改宪法 <- ""
dat$C2a日本修改宪法 <- ""
dat$C2b日本修改宪法 <- ""
dat$C3a日本修改宪法 <- ""
dat$C3b日本修改宪法 <- ""

dat$C1a日本修改宪法[dat$f11=="日本修改宪法"] <- dat$f111[dat$f11=="日本修改宪法"]
dat$C1b日本修改宪法[dat$f11=="日本修改宪法"] <- dat$f121[dat$f11=="日本修改宪法"]
dat$C1a日本修改宪法[dat$f12=="日本修改宪法"] <- dat$f112[dat$f12=="日本修改宪法"]
dat$C1b日本修改宪法[dat$f12=="日本修改宪法"] <- dat$f122[dat$f12=="日本修改宪法"]
dat$C1a日本修改宪法[dat$f13=="日本修改宪法"] <- dat$f113[dat$f13=="日本修改宪法"]
dat$C1b日本修改宪法[dat$f13=="日本修改宪法"] <- dat$f123[dat$f13=="日本修改宪法"]
dat$C1a日本修改宪法[dat$f14=="日本修改宪法"] <- dat$f114[dat$f14=="日本修改宪法"]
dat$C1b日本修改宪法[dat$f14=="日本修改宪法"] <- dat$f124[dat$f14=="日本修改宪法"]
dat$C1a日本修改宪法[dat$f15=="日本修改宪法"] <- dat$f115[dat$f15=="日本修改宪法"]
dat$C1b日本修改宪法[dat$f15=="日本修改宪法"] <- dat$f125[dat$f15=="日本修改宪法"]
dat$C1a日本修改宪法[dat$f16=="日本修改宪法"] <- dat$f116[dat$f16=="日本修改宪法"]
dat$C1b日本修改宪法[dat$f16=="日本修改宪法"] <- dat$f126[dat$f16=="日本修改宪法"]
dat$C1a日本修改宪法[dat$f17=="日本修改宪法"] <- dat$f117[dat$f17=="日本修改宪法"]
dat$C1b日本修改宪法[dat$f17=="日本修改宪法"] <- dat$f127[dat$f17=="日本修改宪法"]

dat$C2a日本修改宪法[dat$f21=="日本修改宪法"] <- dat$f211[dat$f21=="日本修改宪法"]
dat$C2b日本修改宪法[dat$f21=="日本修改宪法"] <- dat$f221[dat$f21=="日本修改宪法"]
dat$C2a日本修改宪法[dat$f22=="日本修改宪法"] <- dat$f212[dat$f22=="日本修改宪法"]
dat$C2b日本修改宪法[dat$f22=="日本修改宪法"] <- dat$f222[dat$f22=="日本修改宪法"]
dat$C2a日本修改宪法[dat$f23=="日本修改宪法"] <- dat$f213[dat$f23=="日本修改宪法"]
dat$C2b日本修改宪法[dat$f23=="日本修改宪法"] <- dat$f223[dat$f23=="日本修改宪法"]
dat$C2a日本修改宪法[dat$f24=="日本修改宪法"] <- dat$f214[dat$f24=="日本修改宪法"]
dat$C2b日本修改宪法[dat$f24=="日本修改宪法"] <- dat$f224[dat$f24=="日本修改宪法"]
dat$C2a日本修改宪法[dat$f25=="日本修改宪法"] <- dat$f215[dat$f25=="日本修改宪法"]
dat$C2b日本修改宪法[dat$f25=="日本修改宪法"] <- dat$f225[dat$f25=="日本修改宪法"]
dat$C2a日本修改宪法[dat$f26=="日本修改宪法"] <- dat$f216[dat$f26=="日本修改宪法"]
dat$C2b日本修改宪法[dat$f26=="日本修改宪法"] <- dat$f226[dat$f26=="日本修改宪法"]
dat$C2a日本修改宪法[dat$f27=="日本修改宪法"] <- dat$f217[dat$f27=="日本修改宪法"]
dat$C2b日本修改宪法[dat$f27=="日本修改宪法"] <- dat$f227[dat$f27=="日本修改宪法"]

dat$C3a日本修改宪法[dat$f31=="日本修改宪法"] <- dat$f311[dat$f31=="日本修改宪法"]
dat$C3b日本修改宪法[dat$f31=="日本修改宪法"] <- dat$f321[dat$f31=="日本修改宪法"]
dat$C3a日本修改宪法[dat$f32=="日本修改宪法"] <- dat$f312[dat$f32=="日本修改宪法"]
dat$C3b日本修改宪法[dat$f32=="日本修改宪法"] <- dat$f322[dat$f32=="日本修改宪法"]
dat$C3a日本修改宪法[dat$f33=="日本修改宪法"] <- dat$f313[dat$f33=="日本修改宪法"]
dat$C3b日本修改宪法[dat$f33=="日本修改宪法"] <- dat$f323[dat$f33=="日本修改宪法"]
dat$C3a日本修改宪法[dat$f34=="日本修改宪法"] <- dat$f314[dat$f34=="日本修改宪法"]
dat$C3b日本修改宪法[dat$f34=="日本修改宪法"] <- dat$f324[dat$f34=="日本修改宪法"]
dat$C3a日本修改宪法[dat$f35=="日本修改宪法"] <- dat$f315[dat$f35=="日本修改宪法"]
dat$C3b日本修改宪法[dat$f35=="日本修改宪法"] <- dat$f325[dat$f35=="日本修改宪法"]
dat$C3a日本修改宪法[dat$f36=="日本修改宪法"] <- dat$f316[dat$f36=="日本修改宪法"]
dat$C3b日本修改宪法[dat$f36=="日本修改宪法"] <- dat$f326[dat$f36=="日本修改宪法"]
dat$C3a日本修改宪法[dat$f37=="日本修改宪法"] <- dat$f317[dat$f37=="日本修改宪法"]
dat$C3b日本修改宪法[dat$f37=="日本修改宪法"] <- dat$f327[dat$f37=="日本修改宪法"]

dat_conjoint$日本修改宪法 <- c(dat$C1a日本修改宪法[dat$f11 != ""], dat$C1b日本修改宪法[dat$f11 != ""],
                        dat$C2a日本修改宪法[dat$f21 != ""], dat$C2b日本修改宪法[dat$f21 != ""],
                        dat$C3a日本修改宪法[dat$f31 != ""], dat$C3b日本修改宪法[dat$f31 != ""])


dat$C1a钓鱼岛主权 <- ""
dat$C1b钓鱼岛主权 <- ""
dat$C2a钓鱼岛主权 <- ""
dat$C2b钓鱼岛主权 <- ""
dat$C3a钓鱼岛主权 <- ""
dat$C3b钓鱼岛主权 <- ""

dat$C1a钓鱼岛主权[dat$f11=="钓鱼岛主权"] <- dat$f111[dat$f11=="钓鱼岛主权"]
dat$C1b钓鱼岛主权[dat$f11=="钓鱼岛主权"] <- dat$f121[dat$f11=="钓鱼岛主权"]
dat$C1a钓鱼岛主权[dat$f12=="钓鱼岛主权"] <- dat$f112[dat$f12=="钓鱼岛主权"]
dat$C1b钓鱼岛主权[dat$f12=="钓鱼岛主权"] <- dat$f122[dat$f12=="钓鱼岛主权"]
dat$C1a钓鱼岛主权[dat$f13=="钓鱼岛主权"] <- dat$f113[dat$f13=="钓鱼岛主权"]
dat$C1b钓鱼岛主权[dat$f13=="钓鱼岛主权"] <- dat$f123[dat$f13=="钓鱼岛主权"]
dat$C1a钓鱼岛主权[dat$f14=="钓鱼岛主权"] <- dat$f114[dat$f14=="钓鱼岛主权"]
dat$C1b钓鱼岛主权[dat$f14=="钓鱼岛主权"] <- dat$f124[dat$f14=="钓鱼岛主权"]
dat$C1a钓鱼岛主权[dat$f15=="钓鱼岛主权"] <- dat$f115[dat$f15=="钓鱼岛主权"]
dat$C1b钓鱼岛主权[dat$f15=="钓鱼岛主权"] <- dat$f125[dat$f15=="钓鱼岛主权"]
dat$C1a钓鱼岛主权[dat$f16=="钓鱼岛主权"] <- dat$f116[dat$f16=="钓鱼岛主权"]
dat$C1b钓鱼岛主权[dat$f16=="钓鱼岛主权"] <- dat$f126[dat$f16=="钓鱼岛主权"]
dat$C1a钓鱼岛主权[dat$f17=="钓鱼岛主权"] <- dat$f117[dat$f17=="钓鱼岛主权"]
dat$C1b钓鱼岛主权[dat$f17=="钓鱼岛主权"] <- dat$f127[dat$f17=="钓鱼岛主权"]

dat$C2a钓鱼岛主权[dat$f21=="钓鱼岛主权"] <- dat$f211[dat$f21=="钓鱼岛主权"]
dat$C2b钓鱼岛主权[dat$f21=="钓鱼岛主权"] <- dat$f221[dat$f21=="钓鱼岛主权"]
dat$C2a钓鱼岛主权[dat$f22=="钓鱼岛主权"] <- dat$f212[dat$f22=="钓鱼岛主权"]
dat$C2b钓鱼岛主权[dat$f22=="钓鱼岛主权"] <- dat$f222[dat$f22=="钓鱼岛主权"]
dat$C2a钓鱼岛主权[dat$f23=="钓鱼岛主权"] <- dat$f213[dat$f23=="钓鱼岛主权"]
dat$C2b钓鱼岛主权[dat$f23=="钓鱼岛主权"] <- dat$f223[dat$f23=="钓鱼岛主权"]
dat$C2a钓鱼岛主权[dat$f24=="钓鱼岛主权"] <- dat$f214[dat$f24=="钓鱼岛主权"]
dat$C2b钓鱼岛主权[dat$f24=="钓鱼岛主权"] <- dat$f224[dat$f24=="钓鱼岛主权"]
dat$C2a钓鱼岛主权[dat$f25=="钓鱼岛主权"] <- dat$f215[dat$f25=="钓鱼岛主权"]
dat$C2b钓鱼岛主权[dat$f25=="钓鱼岛主权"] <- dat$f225[dat$f25=="钓鱼岛主权"]
dat$C2a钓鱼岛主权[dat$f26=="钓鱼岛主权"] <- dat$f216[dat$f26=="钓鱼岛主权"]
dat$C2b钓鱼岛主权[dat$f26=="钓鱼岛主权"] <- dat$f226[dat$f26=="钓鱼岛主权"]
dat$C2a钓鱼岛主权[dat$f27=="钓鱼岛主权"] <- dat$f217[dat$f27=="钓鱼岛主权"]
dat$C2b钓鱼岛主权[dat$f27=="钓鱼岛主权"] <- dat$f227[dat$f27=="钓鱼岛主权"]

dat$C3a钓鱼岛主权[dat$f31=="钓鱼岛主权"] <- dat$f311[dat$f31=="钓鱼岛主权"]
dat$C3b钓鱼岛主权[dat$f31=="钓鱼岛主权"] <- dat$f321[dat$f31=="钓鱼岛主权"]
dat$C3a钓鱼岛主权[dat$f32=="钓鱼岛主权"] <- dat$f312[dat$f32=="钓鱼岛主权"]
dat$C3b钓鱼岛主权[dat$f32=="钓鱼岛主权"] <- dat$f322[dat$f32=="钓鱼岛主权"]
dat$C3a钓鱼岛主权[dat$f33=="钓鱼岛主权"] <- dat$f313[dat$f33=="钓鱼岛主权"]
dat$C3b钓鱼岛主权[dat$f33=="钓鱼岛主权"] <- dat$f323[dat$f33=="钓鱼岛主权"]
dat$C3a钓鱼岛主权[dat$f34=="钓鱼岛主权"] <- dat$f314[dat$f34=="钓鱼岛主权"]
dat$C3b钓鱼岛主权[dat$f34=="钓鱼岛主权"] <- dat$f324[dat$f34=="钓鱼岛主权"]
dat$C3a钓鱼岛主权[dat$f35=="钓鱼岛主权"] <- dat$f315[dat$f35=="钓鱼岛主权"]
dat$C3b钓鱼岛主权[dat$f35=="钓鱼岛主权"] <- dat$f325[dat$f35=="钓鱼岛主权"]
dat$C3a钓鱼岛主权[dat$f36=="钓鱼岛主权"] <- dat$f316[dat$f36=="钓鱼岛主权"]
dat$C3b钓鱼岛主权[dat$f36=="钓鱼岛主权"] <- dat$f326[dat$f36=="钓鱼岛主权"]
dat$C3a钓鱼岛主权[dat$f37=="钓鱼岛主权"] <- dat$f317[dat$f37=="钓鱼岛主权"]
dat$C3b钓鱼岛主权[dat$f37=="钓鱼岛主权"] <- dat$f327[dat$f37=="钓鱼岛主权"]

dat_conjoint$钓鱼岛主权 <- c(dat$C1a钓鱼岛主权[dat$f11 != ""], dat$C1b钓鱼岛主权[dat$f11 != ""],
                        dat$C2a钓鱼岛主权[dat$f21 != ""], dat$C2b钓鱼岛主权[dat$f21 != ""],
                        dat$C3a钓鱼岛主权[dat$f31 != ""], dat$C3b钓鱼岛主权[dat$f31 != ""])



dat$C1a钓鱼岛及周边海域的开发 <- ""
dat$C1b钓鱼岛及周边海域的开发 <- ""
dat$C2a钓鱼岛及周边海域的开发 <- ""
dat$C2b钓鱼岛及周边海域的开发 <- ""
dat$C3a钓鱼岛及周边海域的开发 <- ""
dat$C3b钓鱼岛及周边海域的开发 <- ""

dat$C1a钓鱼岛及周边海域的开发[dat$f11=="钓鱼岛及周边海域的开发"] <- dat$f111[dat$f11=="钓鱼岛及周边海域的开发"]
dat$C1b钓鱼岛及周边海域的开发[dat$f11=="钓鱼岛及周边海域的开发"] <- dat$f121[dat$f11=="钓鱼岛及周边海域的开发"]
dat$C1a钓鱼岛及周边海域的开发[dat$f12=="钓鱼岛及周边海域的开发"] <- dat$f112[dat$f12=="钓鱼岛及周边海域的开发"]
dat$C1b钓鱼岛及周边海域的开发[dat$f12=="钓鱼岛及周边海域的开发"] <- dat$f122[dat$f12=="钓鱼岛及周边海域的开发"]
dat$C1a钓鱼岛及周边海域的开发[dat$f13=="钓鱼岛及周边海域的开发"] <- dat$f113[dat$f13=="钓鱼岛及周边海域的开发"]
dat$C1b钓鱼岛及周边海域的开发[dat$f13=="钓鱼岛及周边海域的开发"] <- dat$f123[dat$f13=="钓鱼岛及周边海域的开发"]
dat$C1a钓鱼岛及周边海域的开发[dat$f14=="钓鱼岛及周边海域的开发"] <- dat$f114[dat$f14=="钓鱼岛及周边海域的开发"]
dat$C1b钓鱼岛及周边海域的开发[dat$f14=="钓鱼岛及周边海域的开发"] <- dat$f124[dat$f14=="钓鱼岛及周边海域的开发"]
dat$C1a钓鱼岛及周边海域的开发[dat$f15=="钓鱼岛及周边海域的开发"] <- dat$f115[dat$f15=="钓鱼岛及周边海域的开发"]
dat$C1b钓鱼岛及周边海域的开发[dat$f15=="钓鱼岛及周边海域的开发"] <- dat$f125[dat$f15=="钓鱼岛及周边海域的开发"]
dat$C1a钓鱼岛及周边海域的开发[dat$f16=="钓鱼岛及周边海域的开发"] <- dat$f116[dat$f16=="钓鱼岛及周边海域的开发"]
dat$C1b钓鱼岛及周边海域的开发[dat$f16=="钓鱼岛及周边海域的开发"] <- dat$f126[dat$f16=="钓鱼岛及周边海域的开发"]
dat$C1a钓鱼岛及周边海域的开发[dat$f17=="钓鱼岛及周边海域的开发"] <- dat$f117[dat$f17=="钓鱼岛及周边海域的开发"]
dat$C1b钓鱼岛及周边海域的开发[dat$f17=="钓鱼岛及周边海域的开发"] <- dat$f127[dat$f17=="钓鱼岛及周边海域的开发"]

dat$C2a钓鱼岛及周边海域的开发[dat$f21=="钓鱼岛及周边海域的开发"] <- dat$f211[dat$f21=="钓鱼岛及周边海域的开发"]
dat$C2b钓鱼岛及周边海域的开发[dat$f21=="钓鱼岛及周边海域的开发"] <- dat$f221[dat$f21=="钓鱼岛及周边海域的开发"]
dat$C2a钓鱼岛及周边海域的开发[dat$f22=="钓鱼岛及周边海域的开发"] <- dat$f212[dat$f22=="钓鱼岛及周边海域的开发"]
dat$C2b钓鱼岛及周边海域的开发[dat$f22=="钓鱼岛及周边海域的开发"] <- dat$f222[dat$f22=="钓鱼岛及周边海域的开发"]
dat$C2a钓鱼岛及周边海域的开发[dat$f23=="钓鱼岛及周边海域的开发"] <- dat$f213[dat$f23=="钓鱼岛及周边海域的开发"]
dat$C2b钓鱼岛及周边海域的开发[dat$f23=="钓鱼岛及周边海域的开发"] <- dat$f223[dat$f23=="钓鱼岛及周边海域的开发"]
dat$C2a钓鱼岛及周边海域的开发[dat$f24=="钓鱼岛及周边海域的开发"] <- dat$f214[dat$f24=="钓鱼岛及周边海域的开发"]
dat$C2b钓鱼岛及周边海域的开发[dat$f24=="钓鱼岛及周边海域的开发"] <- dat$f224[dat$f24=="钓鱼岛及周边海域的开发"]
dat$C2a钓鱼岛及周边海域的开发[dat$f25=="钓鱼岛及周边海域的开发"] <- dat$f215[dat$f25=="钓鱼岛及周边海域的开发"]
dat$C2b钓鱼岛及周边海域的开发[dat$f25=="钓鱼岛及周边海域的开发"] <- dat$f225[dat$f25=="钓鱼岛及周边海域的开发"]
dat$C2a钓鱼岛及周边海域的开发[dat$f26=="钓鱼岛及周边海域的开发"] <- dat$f216[dat$f26=="钓鱼岛及周边海域的开发"]
dat$C2b钓鱼岛及周边海域的开发[dat$f26=="钓鱼岛及周边海域的开发"] <- dat$f226[dat$f26=="钓鱼岛及周边海域的开发"]
dat$C2a钓鱼岛及周边海域的开发[dat$f27=="钓鱼岛及周边海域的开发"] <- dat$f217[dat$f27=="钓鱼岛及周边海域的开发"]
dat$C2b钓鱼岛及周边海域的开发[dat$f27=="钓鱼岛及周边海域的开发"] <- dat$f227[dat$f27=="钓鱼岛及周边海域的开发"]

dat$C3a钓鱼岛及周边海域的开发[dat$f31=="钓鱼岛及周边海域的开发"] <- dat$f311[dat$f31=="钓鱼岛及周边海域的开发"]
dat$C3b钓鱼岛及周边海域的开发[dat$f31=="钓鱼岛及周边海域的开发"] <- dat$f321[dat$f31=="钓鱼岛及周边海域的开发"]
dat$C3a钓鱼岛及周边海域的开发[dat$f32=="钓鱼岛及周边海域的开发"] <- dat$f312[dat$f32=="钓鱼岛及周边海域的开发"]
dat$C3b钓鱼岛及周边海域的开发[dat$f32=="钓鱼岛及周边海域的开发"] <- dat$f322[dat$f32=="钓鱼岛及周边海域的开发"]
dat$C3a钓鱼岛及周边海域的开发[dat$f33=="钓鱼岛及周边海域的开发"] <- dat$f313[dat$f33=="钓鱼岛及周边海域的开发"]
dat$C3b钓鱼岛及周边海域的开发[dat$f33=="钓鱼岛及周边海域的开发"] <- dat$f323[dat$f33=="钓鱼岛及周边海域的开发"]
dat$C3a钓鱼岛及周边海域的开发[dat$f34=="钓鱼岛及周边海域的开发"] <- dat$f314[dat$f34=="钓鱼岛及周边海域的开发"]
dat$C3b钓鱼岛及周边海域的开发[dat$f34=="钓鱼岛及周边海域的开发"] <- dat$f324[dat$f34=="钓鱼岛及周边海域的开发"]
dat$C3a钓鱼岛及周边海域的开发[dat$f35=="钓鱼岛及周边海域的开发"] <- dat$f315[dat$f35=="钓鱼岛及周边海域的开发"]
dat$C3b钓鱼岛及周边海域的开发[dat$f35=="钓鱼岛及周边海域的开发"] <- dat$f325[dat$f35=="钓鱼岛及周边海域的开发"]
dat$C3a钓鱼岛及周边海域的开发[dat$f36=="钓鱼岛及周边海域的开发"] <- dat$f316[dat$f36=="钓鱼岛及周边海域的开发"]
dat$C3b钓鱼岛及周边海域的开发[dat$f36=="钓鱼岛及周边海域的开发"] <- dat$f326[dat$f36=="钓鱼岛及周边海域的开发"]
dat$C3a钓鱼岛及周边海域的开发[dat$f37=="钓鱼岛及周边海域的开发"] <- dat$f317[dat$f37=="钓鱼岛及周边海域的开发"]
dat$C3b钓鱼岛及周边海域的开发[dat$f37=="钓鱼岛及周边海域的开发"] <- dat$f327[dat$f37=="钓鱼岛及周边海域的开发"]

dat_conjoint$钓鱼岛及周边海域的开发 <- c(dat$C1a钓鱼岛及周边海域的开发[dat$f11 != ""], dat$C1b钓鱼岛及周边海域的开发[dat$f11 != ""],
                        dat$C2a钓鱼岛及周边海域的开发[dat$f21 != ""], dat$C2b钓鱼岛及周边海域的开发[dat$f21 != ""],
                        dat$C3a钓鱼岛及周边海域的开发[dat$f31 != ""], dat$C3b钓鱼岛及周边海域的开发[dat$f31 != ""])

dat_conjoint_sub <- dat_conjoint[dat_conjoint$id %in% dat$id[which(dat$q31=="其他")], ]
dat_conjoint_covariate <- dat[, c(1:which(colnames(dat) == "q310"), which(colnames(dat) == "id"):which(colnames(dat) == "fight"))]

save(dat_vignette, dat_vignette_sub,
     dat_conjoint, dat_conjoint_sub, dat_conjoint_covariate, 
     file="data/dat_analysis_CN.RData")

